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ABSTRACT 

cz) ! 

^ , Context. The determination of fundamental parameters for horizontal hranch stars in Galactic globular clusters frequently gives puzzling re- 
sults, pointing to the lack of our understanding of their atmospheric structure and the inadequate approximations by the models. 
£SJ ■ Aims. We wanted to measure effective temperatures, surface gravities, helium abundances, and masses for hot horizontal branch stars in 

^ ' NGC6752 in order to compare the results with evolutionary predictions. 
t-H Methods. We used the ESO VLT-FORS2 facility to collect low-resolution spectra of 51 targets distributed along the horizontal branch. We 
determined atmospheric parameters, by comparison with theoretical models through standard fitting routines, and masses by basic equations. 
Results. Results generally agree with previous works, although not always with the theoretical expectations for cooler stars (T c ff < 15 000 K). 
' I , The calculated color excess is systematically lower than the literature values, pointing towards a possible underestimation of effective temper- 
OO 1 atures. Moreover, we find two groups of stars at T e j ~ 14 000 K and at T e ff ~ 27 000 K that present anomalies with respect to the general trend 
CZ^ and expectations. We suppose that the three peculiar bright stars at T e j ~ 14 000 K are probably affected by an enhanced stellar wind. For the 
C i eight extreme horizontal branch stars at T c ff ~ 27 000 K that show unusually high masses, we find no plausible explanation. 

' Conclusions. While most of our results agree well with the predictions of standard horizontal branch evolution, we still have problems with the 
^ low masses we derive in certain temperature ranges. We believe that Kurucz ATLAS9 LTE model atmospheres with solar-scaled abundances 
• »"H . are probably inadequate for these temperature ranges. Concerning the group of anomalous stars at T cff ~ 27 000 K, a Kolmogorov-Smirnov 
' test indicates that there is only an 8.4% probability that these stars are randomly drawn from the general distribution in the color-magnitude 
diagram. This is not conclusive but indicates that these stars could be both (and independently) spectroscopically and photometrically peculiar 
with respect to the general extreme horizontal branch population. 
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1. Introduction 

Many decades have passed since when horizontal branch 
(HB) stars were identified as evolved stars of low initial 
mass (app roximately 0.7-2 M^) presently burning heli um in 
their core dHoyle & Schwarzschildl955l:lFaulknerll966l) . Still, 
many aspects of the evolution and internal structure of these 
stars are not fully understood. In particular, in the color- 
magnitude diagrams of Galactic globular clusters, many well- 
established observed features of HB stars are good example 
of this, as seen in the case of the different blue extensions 
from one clus ter to another, which are p artially dependent 
on metallicity (San dage & Waller stein 1960) but not fully ex- 



* Based on observations with the ESO Very Large Telescope at 
Paranal Observatory, Chile (proposal ID 69.D-0682) 



plained by it dSandage & Wildevll 19671 Ivan den Berghlll9671) . 
In addition, the HB may show "jumps ", where stars appear 
bright er than the theo r etical expectations dGrundahl et al . 1999 
iMomanv et al.l 12002, 200 4j), or gaps, i.e . underpopula t ed re 



gions along the HB Tsosin et all [l 9971; ICatelan et all Il998 



Ferraro et al.1 1 19981: IPiotto et al.l Il999l) . Spectroscopy of HB 



stars opened new questions, such as the unexplaine d presence 
of fa st rotator s among stars redder t han th e Grundahl et al 
d 1999b iump JPetersorjl 1 98311 1 985alTbi: lBehrll2003l: iBehr et al 



2000ab ; Recio-B lanco et al.l 120044) and the absence of close 
binaries among hot extr eme horizontal branch (EHB) stars 
dMoni Bidin et al.ll2006alfbh . 

Traditionally, low-resolution spectroscopy of HB stars has 
focused on determining of their fundamental quantities such 
as atmospheric parameters. This is usually achieved by com- 



2 



Moni Bidin et al.: Spectroscopy of horizontal branch stars in NGC6752 



paring observed spectra with grids of theoretical ones obtained 
from model atmospheres. As a consequence, any aspect of stel- 
lar structure unaccounted for in the model may be directly re- 
flected on unexpected values measured for the parameters. For 
example, it is well known that surface gravities are systemat- 
ically lower than expected in the temperature range 12000- 
20000 K if model spectra with the metallicity of the parent 
globular cluster are us ed. As a consequenc e, masses are de- 



rived that are too low (Moehl er et al. I ll995l) . This "mass dis 



crepancy" can b e partially explained b y radiative levitation of 
heavy elements ( Grundahl et al J [19991) . originally proposed as 
the physical process res posible for ab undance anomalies ob- 
served in field Ap stars dMichaudlll970l) . Precise calculations 
confirmed that this phenomenon should be at work in the at- 
mos pheres of these stars a nd can account for observed anoma- 
lies (Mi chaud et al.lll983l) . but the pr oblem of the m ass dis- 
crepancy is still not completely so lved ( Moehlerll2001 ). On the 
other hand. lVink & Cassisil ( 2002) point out that neglecting the 
presence of stellar wind can cause measured surface gravities 
to be erroneously low. Nevertheless, the derived values of stel- 
lar parameters - even if erroneous or unexpected - have been 
useful for pointing out where canonical models are inadequate, 
stimulating further investigations to solve the open problems. 

As part of a pr ogram aimed at searchin g for close binaries 
among HB stars dMoni Bidin et al.ll2006bl hereafter Paper I), 
we also collected low-resolution spectra of our target stars. In 
this paper we present the analysis of these spectra and critically 
discuss the derived atmospheric parameters (effective tempera- 
ture, surface gravity, helium abundance) and stellar masses. 



2. Observations and data reduction 

The spectra were collected during two nights of observations 
(June 12 and 13, 2002) at the VLT-UT4 telescope equipped 
with the spectrograph FORS2 in MXU mode. The fifty-one 
target stars are the same objects already described in Paper 
I. They are distributed along the entire HB, spanning a wide 
range in temperature from T e ff ~ 8 QQQK up to T e ff ~ 30 QQQK. 
Coordinates and photometric data of the observed stars are pre- 
sented in Table 1 of Paper I. Figure Q] shows the position of 
target stars in the color-magnitude diagram. 

Two 1350s exposures were secured for each target with 
grism 600B+22 and 075 - wide slits for a resulting resolution 
of 3 A. The spectral range was approximately 2900 A wide 
and, on average, centered at 4600 A, but slightly different from 
star to star due to the different positions of the slits in the mask 
used for multi-object spectroscopy. All Balmer lines from Hg to 
H10 were always present in the spectra. We never used wave- 
lengths shorter than 3600 Ain our analysis, because of the lack 
of atmospheric transmission and instrumental response. In Fig. 
|2]we show an example of the reduced spectra. 

The calibration images (bias, flat, and lamp) were acquired 
before and at the end of each observing night. Data reduction 
was performed with standard MIDASLj procedures. All of the 
2D spectra were trimmed from the frames with their calibra- 
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Fig. 1. Observed stars in the V vs U — V color- magnitu de di- 
agram. Photometric data are from iMomanv et aL ( 2002 ). The 
stars are indicated with different symbols, as discussed in the 
text. The same symbols will be used in all of the figures to dis- 
tinguish groups of stars with different properties. 



tion images and reduced independently. The wavelength cal- 
ibration was performed with the HeHgCd lamp images, fit- 
ting a 3 rd order polynomial to the dispersion relation. The rms 
error indicated by the fitting procedure was 0.07 Aon aver- 
age. Before extraction, we rebinned the 2D images to constant 
wavelength steps of 0.4 A/pix and corrected their curvature 
along the spatial axis by t racing the spectra as described in 
iMoehler & Sweigarj (120061) . Spectra were then ext racted both 
with an optimum extraction algorithm ( Hornd 19 86) and with a 
simple sum, and then, for each star, we chose the procedure that 
gave better results in terms of S/N, as the optimum algorithm 
sometimes failed to correctly extract the spectra of bright stars. 
We corrected them for atmospheric extinctio n, with the extin- 
tion coefficients for the La Silla observatory (Tug 1977). The 
response curve was obtained with observations of the standard 
stars EG274 (both nigh ts) and LTT3218 (s econd night only), 
using the flux tables of Ham uv et al. (1994). We had to use a 
high-order polynomial (6th) so as to obtain a good response 
curve and a smooth continuum on the calibrated spectra. The 
resulting curves had the same shape with no difference between 
the nights. To calibrate the spectra of the second night we av- 
eraged the two corresponding curves when two standard stars 
were observed. 

Finally, we corrected the spectra for radial velocity shifts. 
Absolute radial velocities of the targets were already measured 
on higher resolution spectra as described in Paper I, but we 



1 ESO-MIDAS is the acronym for the European Southern 
Observatory Munich Image Data Analysis System which is de- 



veloped and maintained by the European Southern Observatory 
( http://www.eso.org/projects/esomidas/ 1 
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Fig. 2. A sample of the collected spectra, normalized and ordered by increasing T e ff from top to bottom. Spectra were normalized 
fitting a 3 rd - order spline function to the continuum. For wavelengths shorter than 3900 A, where Balmer lines blend, the fit was 
anchored at the highest points between lines. 



measured them again on each of the low-resolution spectra 
used in this work. In fact, as already discussed in Paper I, the 
shift of the spectra with respect to laboratory wavelength was 
due not only to the radial velocities of the stars, but also to 
many effects such as the position of the stars in the slitlets and 
systematics in wavelength calibration. We fitted the cores of all 
Balmer lines from H^ to H9 with a Gaussian profile, exclud- 
ing H e due to blending with the CanH line, and assumed the 
average as the radial velocity to be corrected. 

Usually the values from single lines agreed within 10 km 
s _1 , but in some cases there were large diffe rences (up to 5 km 
s~'), an effect already described in iMoehler et alJ (11997b . We 
checked if these could be ascribed at least in part to a rigid shift 
of the spectra on the CCD, constant in pixel and then variable 
with wavelength in km s , but we found no evidence of a clear 
trend with wavelength. 



3. Measurements 

In order to derive effective temperatures, surface gravities, and 
helium abundances, we fitted the observed Balmer and Helium 
lines with stellar model atmospheres. 

As discussed extensively in recent years, HB stars hot- 
ter than about 11 000 - 12 000 K show evidence of diffusion, 
e.g. deficiency of helium and strong (solar to super-solar lev- 
els) e nrichment in iron (fGrundahl et alJ Jl999t IGlaspev et al, 
1 19851 NGC 6397: IGlaspev et alJ Il989r. IMoehler et all l200l 



NGC6752 ; iBehr etalJ 1 19991: IMoehler et alJ 120031 M13; 



Behretal. 2000a, M15; Fabb ian et all 120051 NGC1904; 



v- 
tl. 
0, 



Pace et alJ 120061 NGC2808). As the atmospheric metallicity 



49317 (7800 K) 



4300 




4600 



4700 



has a non-negligible influence on the profiles of hydrogen and 



Fig. 3. A sample of spectra for stars across the Grundhal jump 
in the spectral range used to determine the metallicity of the 
model in the fitting procedure. Some of the most relevant Fe n 
are indicated with the Mg 11 and He 1 doublets that are the most 
important features in this spectral range. 



helium lines, we need to know if diffusion is active in a star 
to decide on the metallicity to fit our spectra. Therefore, in all 
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spectra we carefully examined the region 4450 A to 4600 A, 
where strong Fe n lines show up at 1 1 00 K to 12 000 K. 
Past experience dMoehler et al.ll 19991 12000[ 120031) has shown 
that these lines are strong enough to allow detection even 
in medium-resolution spectra. In Fig. [3] we show this wave- 
length range for some spectra of stars across the Grundhal 
jump. Spectra showing evidence of iron lines or of being hot- 
ter than 14000 K (deduced from their position in the color- 
magnitude diagram) were fit with metal-rich ([M/H] +0.5) 
model spectra, whereas stars below 14000 K with no evidence 
of strong iron lines were fit with metal-poor model [M/H] = 
—1.5) spectra. In the cool star spectra, we kept the helium 
abundance fixed, as the helium lines in cool stars are rather 
weak, and the helium abundance should be close to solar any- 
way. During the fitting we verified that those helium lines pre- 
dicted for these cool stars agree with the obse rved ones. We 
computed model atmospheres using ATLAS9 dKuruczlll993l) 
and used Lemke's versiorQ of the LINFOR program (devel- 
oped originally by Holweger, Steffen, and Steenbock at Kiel 
University) to compute a grid of theoretical spectra that include 
the Balmer lines H a to H22, He 1, and He n lines. The grid cov- 
ered the range 7000 K < T eff < 35 000 K, 2.5 < logg < 6.0, 
-3.0 < log He < -1.0, at metallicities of [M/H] = -1.5 and 
+0.5. In Table Q] we list the results obtained from fitting the 
Balmer lines Hp to H10 (excluding H e to avoid the Can H line) 
and, in the hot stars, also the He 1 lines 4026 A, 4388 A, 4471 A, 
4921 A and the Hen lines 4542 A and 4686 A. IDs are from 



Table 1. Atmospheric parameters and derived masses for target 
stars. 



Mo many et alJ (120021) . In cool star spectra the He abundance 



was kept fixed (-1.00+0.00). Errors comes from the weighting 
procedure of two measurement for each star, and are multiplied 
by a/3 but not corrected for possible underestimation. Masses 
were derived with Eq. Q] In last column peculiar stars are indi- 
cated, with the same symbols used in all figures. 

To establish the bes t fit to the observed spe ctra, we used the 
routin es developed by iBergeron et alJ dl992h and Saffer et al 



( 1994), as modified by Napi wotzki et al.l (119991) . which employ 
ax 2 test. The <x necessary for the calculation ofx 2 is estimated 
from the noise in the continuum regions of the spectra. The fit 
program normalizes model spectra and observed spectra using 
the same points for the continuum definition. 

For each star, the atmospheric parameters were measured 
independently in the two low-resolution spectra collected dur- 
ing the observing run. The final results given in Table Q] are 
the weighted means of these values, weighted by the inverse 
errors provided by the fit procedure. The tabulated errors are 
the results of the weighted mean procedure applied to the pair 
of independent measurements obtained by the fit of the two 
spectra for each star. We multiplied the errors given by the fit- 
ting routine by V3 because, when rebinning, we oversampled 
the spectra by a factor of three with respect to the dispersion, 
while the fitting procedure assumes each pixel is independent 
of the others. The errors in each fit of the observed spectrum 
with the model spectrum were derived from the x 2 of the fit 
itself (see Mo ehler et al. 1999, for more details). These errors 
are obtained under the assumption that the only error source 
is the statistical noise (as derived from the continuum of the 
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spectrum). However, Napiwotzki (priv. comm.) noted that the 
routine underestimates this statistical error by a factor of 2-4. 



For a description see http://a400.sternwarte.uni-erlangen.de/~ai26/linfit/linfor.html 
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In addition, errors in the normalization of the spectrum, im- 
perfections of flat field/sky background correction, etc. may 
produce systematic rather than random errors, which are not 
represented well by the error obtained from the fit routine. 
Masses were calculated from the previously measured atmo- 
spheric parameters in two ways. First we used the equation 

M 

log = const. + \ogg + 0.4 • ((;« - M) v -V- V th ) (1) 

M Q 

where V t h is t he brightness at the stellar surface as given by 
Kuruczl d 1992b . In addition we also calculated masses through 
the equation 



Ms T 
log — = log 4 • log — 

Mq go Tq 



log 



Lq 



(2) 



obtained from basic relations, assuming Ilq-5111 K and 
logg0=4.4377. In the last term we expressed the luminosity 
in terms of absolute magnitude and bolometric correction BC 
through the relation 



log- 



My + BC- 4.74 
Z5 ' 



(3) 



The bolometric correction for each target star was derived 
from ef f ective temperatures using the empirical calibration of 



Flower ] (| 1996b. In both eq uations we used photometric data 
from llVIomany et aL ( 20021) . and we adopted an apparent dis- 



tance modulus (m-M)y=13. 17 and Eb-v=0.04 as the mean 
value of the determination of Renzini et al. d 1 9961) . iReidl (1 19971 
1998b and iGratton et ail dl997t) . Results from both equations 
are plotted separately in Fig. [7] Errors on masses were derived 
from propagation of errors, including the V3 factor for errors 
on T e ff and logg. We also assumed an error of 0. 1 magnitudes 
on the photometric quantities used in the equations, such as 
distance modulus, magnitude, V,/,, and bolometric correction. 

4. Results 

4.1. Temperatures and gravities 

Our results on atmospheric parameters are plotted in Fig. H] 
where the position of the target stars in the T e ff-logg space is 
compared with theoretical models. The agreement between the 
measured parameters and the expectaction from the models is 
good, and all stars but a few exceptions fall between the theo- 
retical zero-age HB (ZAHB) and the terminal-age HB (TAHB) 
for normal He-content (Y=0.24) stars. 

We note a group of hot stars (T e ff ~27 000 K) for which the 
measured surface gravities seem to be systematically higher 
than expected for normal He content stars. Another group of 
three cooler stars (T e fr~14000 K) shows gravities that are too 
low and lies above the TAHB. Though these differences are 
marginal in the T^-logg plot, they appear more significant 
when we consider the masses (see Fig. [7J. It is worth notic- 
ing in Fig |4] that they are all reproduced better by models with 
enhanced helium abundance (Y=0.33). 

There is an additional anomalous star at T e ffa;25 000 K, 
which has a significantly smaller surface gravity than expected. 
This star (marked as asterisk in all the figures in this paper) also 
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Fig. 4. Temperatures and gravities of our program stars. All 
stars hotter than 12 000 K or showing evidence of iron enrich- 
ment are fitted with model atmospheres of super-solar metal- 
licity ([M/H] = +0.5). Cooler stars showing no evidence of 
iron enrichment are fitted with metal-poor model atmospheres 
([M/H] = - 1 .5). The solid lines mark the zero-age (ZAHB) and 
terminal-age (TAH B) loci of canonical H B tracks (Y=0.24) for 
[M/H] = - 1 .54 (see lMoehler et al.l2003l for details). These loci 
define the region within which the HB models spend 99% of 
their HB lifetime. Dashed lines represent the ZAHB and the 
TAHB of helium enriched (Y=0.33) stars. Error bars are omit- 
ted for clarity, but the typical size of the errors as derived from 
the x 1 fit is indicated. The thick error bars indicate the formal 
error as in Table [T] and the larger thin error bars indicate three 
times the previous value for hotter stars, approximately 1000 K 
and 0. 1 dex, which is probably a better estimate of real errors 
(see text for details). 



has an anomalous position in the color-magnitude diagram and 
will be discussed later. All the anomalous stars have radial ve- 
locities compatible with the radial velocity of NGC6752, as 
can be seen in Table 1 of Paper I. Therefore they can be con- 
sidered probable cluster members, but it must be noted that ra- 
dial velocity alone is not an unambiguous proof of membership 
in this cluster. In fact, at the galactic coordinates of NGC6752 
(Z=336°, b- -26°), a certain amount of field contamination is 
expected, and its radial velo city is unrem arkable compared to 
field disk stars (-27.9 km s~'. lHarrisll996h . 

We comp ared the derive d effe ctive temperatures with those 
measured in Moe hler et al. ( 2000h for HB stars in the same 
cluster. The temperatures of both data sets are plotted as a func- 
tion of color in Fig. [5] The results agree very well. There is a 
small offset of our data toward lower temperatures, but it cannot 
be considered significant because it is much smaller (< 500 K) 
than formal errors. We also have nine stars in common with 
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Fig. 6. Differences i n surface gravities and derived masses be- 
tween this work and Moehler et al. (2000) for the nine stars in 
common. The errors are the quadratic sum of the ones for each 
set of data. The difference is in the sense (ours) - (MOO). 



Fig. 5. Measured effective temperatures as a function of colors 
derived in this work (full circles, and triang les and squares fol- 
low the notation in previous plots) and in lMoehler et al.l J2 000) 
(small empty circles). The adopted data from iMoehler et al.l 
(2000) are the ones obtained with metal-poor model atmo- 
spheres (their Table 2) for stars with T e jf < 1 1 000 K and with 
metal-rich model atmospheres (their Table 5) for hotter stars. 



Moe hler et al. ( 2000l) . In Fig.[5]we plot the differences between 
the two works in derived logg and masses. Within the errors, 
the results of the two investigations are consistent, as was ex- 
pected due to the similarities between the two works: data re- 
ductions are almost identical, and measurement procedures and 
model atmospheres are the same. 

4.2. Masses 

In Fig. [7] we plot the masses we derived for our target stars as 
a function of effective temperatures, with both Eqs. [TJ and [2] 
Here the problems and the differences between observed and 
expected values become evident. We note that both equations 
give very similar results, and the differences are negligible in 
all but the temperature range 12000-15 000 K. In this interval 
the masses calculated with Eq. [2] show better agreement with 
theoretical expectations because they are on average 0.05 M0 
higher, but it is not enough to correct the general tendency of 
an underestimation of the derived masses for T e ff < 15 000 K. 
We discuss this problem in more details in subsequent sections. 

Here we note that, although the plotted points in this tem- 
perature range are closer to the theoretical line in the lower 
panel of Fig. [7] they are not scattered around it, but system- 
atically below it, as in the upper panel. This indicates that the 
problems concerning derived masses cannot be ascribed only 
to either the theoretical V,/, values used in Eq. [TJor to empiri- 



cal BCs used in Eq.[2] because both equations lead to the same 
underestimation independently. 

The problem cann ot be attributed to the V magnitudes from 
Momanv et alj (12002b that we use to derive the masses. When 
calibrating photometric data, we all refer to the standard UBV 
Johnson filters and calibrate w ith respect to them. However, the 
study bv lMomanv et al.l (2003) basically shows that no 2 filter- 
systems are the same (see their Fig. 2) and that small differ- 
ences in the bandpass of the employed filters (at any telescope) 
may persist. Momany et al. also show that such "small" dif- 
ferences can eventually give rise to un-physical features, and 
we are aware that the WFI@2.2m V filter has a red cutoff at 
A ~ 5900A with respect to the standard V filter. Moreover, the 
WFI@2.2m instrumental v magnitudes were calibrated with 
their respective (b— v) colors, thus employing the WFI@2.2m B 
filter with a redder effective central wavelength. Nevertheless, 
it remains quite hard to explain the derived low masses in 
terms of photometric calibration uncertainties. Indeed, to de- 
rive masses that match theoretical expectations, we estimate 
that the My magnitudes need to be increased by ~ 0.4 in the 
12 000 - 15 000 K range, and decreased by a similar amount 
for the hot outliers. Clearly, this discrepancy is too large to be 
attributed to photometric calibration uncertainties. 

Because of the lack of noticeable differences between the 
panels in Fig. [7] we discuss the results without distinguishing 
between the two methods to derive masses. Instead, we find 
it appropriate to divide the sample to three different ranges of 
temperatures. 

4.2.1. Masses: cool stars (T eff <12 000 K) 

The derived masses of stars cooler than 10000 K are sys- 
tematically too low. All of these st ars are offset from the HB 
towards lower gravities in Fig. [4] Moehler et al. (2003) ob- 
served the same effect for HB stars in Ml 3 and excluded these 



stars fr om their analysis. It is worth noticing that Moeh ler et al 



(2000) find masses in this range of temperatures in good agree- 
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ment with the expected valu e for stars in the same cluster, but 
Moehler & Sweigartl d200fj|) in NGC6388 find masses that are 



too low, as we do with the same instruments and set-up as in 
the present work (UT4+FORS2 and grism B600). They con- 
sider that atmospheric parameters for these stars are probably 
not trustworthy, because of the presence of metallic lines not 
considered in the fitting procedure and not clearly observable 
due to low resolution. In the case of NGC6752, metal lines offer 
less explanation because this cluster is more than a factor of 10 
more metal-poor than NGC6388. However because of the co- 
incidence observed, it can be hypothesized that, whatever the 
cause, it could somehow be enhance d by the instrumentation 
and set-up used. lMoehler et al.1 (120031) find too low mass for the 
star they observed in this temperature range, in a cluster of sim- 
ilar metallicity (M3) but with different instrumentation. This 
excludes the possibility that the underestimation was only due 
to instrumental effect with this set-up, and the problem of the 
masses too low for T e ff < 10 000 K remains completely open. 
The calculated masses for stars between 10000 and 12 000 K 
better agree with theoretical values although the scatter is quite 
high. 



4.2.2. Masses: intermediate temperature stars 
(12 000<T eff <15 000 K) 

For higher temperatures, as the radiative levitation sets in, 
masses are constantly underestimated. This is a well-known re- 
sult that is clearly present in our data. In fact, even if the errors 
on the masses of the single stars are high, Fig. [7] clearly shows 
that our derived masses in the 12 000<T e ff <15 000 K interval 
are systematically lower tha n predicted. Our resu lts are also in 
agreement with the ones of iMoehler et alj (2000), and the dif- 
ference that we find between observed and expected masses is 
of the same order of magnitude as in their work. We then con- 
firm their conclusions that the use of models with higher metal- 
licity in the fitting procedure, to account for radiative levitation, 
can only partially reduce the discrepancy. 

There is a group of three stars located at 13 000 <T e ff < 
14 000 K with particularly low masses and low gravities. These 
stars are indicated in all of the figures as empty triangles. These 
stars are also separated well in the color-magnitude diagram 
(Fig. Q]), as they are systematically brighter than the others. It 
must be emphasized that, although we relied on photometric 
data to derive masses, the difference in masses cannot be as- 
cribed to the difference in magnitude, since they go in oppo- 
site directions: all other parameters being the same, brighter 
stars should show higher masses, and not lower, as can be de- 
duced by Eq. Q] These anomalous masses are probably due to 
the very low surface gravities, low enough to give underesti- 
mated masses, even if the stars are brighter than their neighbors 
in the color-magnitude diagram. Therefore, for some reason, 
these stars have different luminosities and a different spectrum 
than other HB stars of similar temperatures. In this temperature 
range, there is another star (#18782) brighter than the average 
HB in Fig. [T]and outside (above) the theoretical tracks in Fig. 
|U Nevertheless its calculated mass is not noticeably low with 
respect to other stars with similar temperatures, so we did not 
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Fig. 7. Estimated masses of target stars as a function of effec- 
tive temperature. Different symbols are used to distinguish be- 
tween groups of stars that behav e differently, as ind icated in 
the text. Also the theoretical HB (M oehler et al]l2003l) is indi- 
cated. Error bars on temperatures are omitted for clarity, but are 
within an order of magnitude of full points. Errors on masses 
are obtained from the propagation of errors on temperature and 
gravities, this error bars referring to the formal errors given in 
Table Q] and larger ones assuming errors of 1000 K and 0.1 
dex respectively, which are probably more realistic estimates 
of real errors. Vertical dashed lines divide the temperature in- 
tervals analyzed in the text. Upper panel: masses derived from 
Eq.Q] Lower panel: masses derived from Eq.|2] 



point it out with a different symbol. We do take it into account 
in the later discussion. 



4.2.3. Masses: hot stars (T eff >15 000 K) 

For increasing temperatures, the measured masses again 
generally agree with theoretical values, as also found by 
Moehler et alj (120001) in the same cluster. Nevertheless, also 
among hotter stars there is a group of stars that behave dif- 
ferently from the general trend. These objects are indicated as 
empty squares in all the plots, and for them the derived masses 
are systematically higher than expected. As shown in Fig. Q] 
they are fainter than the bulk of HB stars at the same tempera- 
ture. Again, Eq.Q] would imply lower masses for fainter stars, 
so this is not caused by how we derived our masses: either these 
objects are intrinsically more massive or they are both photo- 
metrically and spectroscopically different from other stars at 
the same temperature. 
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Fig. 8. Derived Helium abundances as a function of effective 
temperatures for target stars. Different symbols are used to dis- 
tinguish between groups of stars that behave differently, as in- 
dicated in the text. The dotted line indicates solar value, and the 
typical errorbars (0.09 dex in log(N((He)/N(H)) and 300 K in 
temperature) are indicated in upper right. For stars cooler than 
the Grundhal jump the abundance was kept fixed to solar value 
in the fitting procedure, and they are omitted in the plot. 



4.3. Helium abundances 

The derived He abundances are plotted in Fig.|8]as a function of 
effective temperatures, where only stars hotter than Grundhal 
jump are considered because, for cooler ones, He abundance 
was not measurable, hence kept fixed to solar value. All stars 
show helium depletion to a different extent, but no trend with 
temperature is visible. 

The three stars indicated as triangles, which show particu- 
larly low gravities and/or masses, on average also show a higher 
helium abundance than do other stars at the same temperature, 
though they agree with values for stars at different temperatures 
- from cooler ones up to 20 000 K. On the other hand, no clear 
trend is visible for hot stars with anomalous measured masses, 
which show helium abundances in agreement with other hot 
stars. 

4.4. Star #17941 

Star #17941 is very peculiar and its measured atmospheric pa- 
rameters are highly uncertain. This object is indicated with an 
asterisk in all the figures. 

Its radial velocity, as measured in Paper I, indicates that it is 
a probable cluster member. In the color-magnitude diagram it is 
located below the HB at fainter magnitudes and/or redder color. 
Its B-V color indicates an effective temperature around 14000 
K, which is also confirmed by its intermediate-resolution spec- 
trum (Paper I), which shows the typical features of an HB star 
with this temperature, in particular, a high quantity of metallic 
lines due to radiative levitation. Actually, its spectrum is very 
similar to the ones of other stars in this temperature range. On 
the contrary, from the measurements on low-resolution spectra, 
it comes out to be a much hotter HB star (24800 K). Also the 
surface gravity is peculiar, being particularly low. Even more 
puzzling, and at variance with the behavior of all other pecu- 
liar stars already analyzed, we find a higher mass than expected 
even if its surface gravity is lower than canonical models. 
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Fig. 9. Estimated E(b-v) reddening values for target stars. The 
dotted line indicates the valu e E(b-vo=0-04, whic h is the mean 
value of the determination of Renzini et al. dl996l) . lReidldl997l 
1998), and lGratton etail dl997l) . 



We analyzed the data reduction of these spectra in de- 
tails, trying different solutions for spectra reduction and ex- 
traction, with no significant changes in the derived parameters. 
We found no evidence of light contamination by a nearby star, 
both when analyzing the slit images and in the field images ac- 
quired before exposures. We did not find any indication of any 
companion, not even among the spectral features observed at 
intermediate resolution, although this could not be indicative 
because the resolution itself is not very high (1.2 A), and the 
great quantity of metallic lines could easily hide fainter fea- 
tures from a companion. 

4.5. Reddening 

With the derived effective temperatures and gravities, we es- 
timated the color excess E(B-V) that could be deduced from 
these atmosph eric parameters. We used the Kurucz ATLAS9 
dKuruczll 993) photometry table for the same metallicity as was 
used in fitting procedure, and for each star we interpolated the 
B-V values in the grid to get the expected values. Then we cal- 
culated the color eccess E(B-V) b y comparing this theor etical 
color with the observed one from iMomany et aL (2002) pho- 
tometry. The results are plotted in Fig. [9] 

The computed reddening is constantl y too low compared t o 
litera t ure values. I n fact the mean value of | Renzini et alJ ( 19961) . 
Reidl d 19971 1 1998b . and lGratton eTail dl997b that we also used 
in Eq. Q]is E(B-V)=0.04, an d it is in fair agreem ent with the 
more recent measurement bv lGratton et alj (120051) . We obtain 
E(B-V)=0.016 as the mean value for our stars, with no partic- 
ular trend with temperature. 

It should be noted that the WFI@2.2m bandpasses dif- 
fer significantly from the standard Johnson-Cousins ones 
dMomany et a il l2002h . and an error of 0.02 magnitudes in 
color cannot be excluded. However. iGratton et al.l (12005) use 
the same photometric data and obtain reddening estimates in 
perfect agreement with previous works among turn-off stars 
in the same cluster. The main difference is that they esti- 
mated temperatures fitting H ff profiles in high-resolution spec- 
tra, a totally different procedure from the one adopted here. 
Therefore, we are forced to conclude that the offset observed 
here could probably point to a systematic underestimate of 
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the temperature along the HB. Of course, the problem could 
reside in the theoretical colors themselves. In fact, empirical 
calibrations of T e ff-color relations show that synthetic colors 
from model atmospheres calculated wit h ATLAS9 have diffi- 
culty properly reproducing observat ions dRamrrez & Melendezl 
20051: ISekiguchi & Fukugital Eoool) . a proble m also pointe d 
out in the theoretical investigation by ICastelli et alJ dl997l) . 
Unfortunately, all the cited empirical calibrations are limited 
in the cool temperature range (spectral type G-K) of main 
sequence or giant stars, and there is a general lack of HB 
color-temperature, model-independent calibrations in the liter- 
ature. In conclusion, it is not possible to guess to what extent 
ATLAS9 models could lead to a systematic error in tempera- 
ture. 

The stars that show peculiar masses show no particular be- 
havior in Fig. [9] Since reddening estimates are very sensitive 
to temperature but should hardly change for small variations of 
surface gravity, this can indicate that their derived temperature 
are not peculiar, and their strange behavior is only due to an un- 
usual surface gravity (either real or as a consequence of other 
effects unaccounted for). 

5. Discussion 

We assumed LTE in our model spectra, as usual in this kind of 
analysis, but it is becoming evident that non-LTE effects could 
affect the results for hotter stars and lead to erroneous a tmo- 
spheric parameters, as pointed out by lPrzybilla et al. ( 2006h . In 
their preliminary results they show that a hybrid NLTE analysis 
can provide a much better fit of H and He lines, but the differ- 
ences in atmospheric parameters relative to results assuming 
LTE are only slightly larger than our real uncertainties (on the 
order of 1000 K in T e ff and 0.1 dex in logg). Considering that 
errors in Table Q] are underestimated by a factor of at least 2-4, 
as already discussed, the influence of neglecting NLTE effects 
is relatively small. Even more important, it should affect all 
stars in the same way at a given temperature (NLTE analysis 
provides higher temperatures and gravities than the LTE one), 
so it should not cause any systematic difference for stars at ap- 
proximately the same temperature. However, adding 1000 K in 
T e ff and 0. 1 dex in logg strongly increases the derived masses 
as calculated from Eq.Q] The change in temperature is reflected 
in a decrease of ^0.1 mag for V,/,, which is quite insensitive to 
gravity. Masses are then higher by about 0.2 Mq in the high 
temperature range. The effect on masses derived from Eq.|2]is 
much smaller, approximately half the previous one, but again in 
the sense of increasing masses. While this would alleviate our 
problems between 12 000 K and 15 000 K, there is no reason to 
assume that NLTE effects would be limited to this temperature 
range. 

We also assumed solar abundance ratios in the models. 
This is necessarily a crude approximation in order to take into 
account the enhanced abundancies of heavy elements due to 
radiative levitation. Detailed st udies both in the p roximity of 
the Grundhal Jum p dBehr et al.ll 19991: lBehrll2003l) and among 
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hotter HB stars ( Edelma nn et al J 1200 1 ) show that the atmo 



spheric metal abundances in the presence of radiative levita- 
tion do not usually follow solar ratios: the elements are af- 



Fig. 10. A section of the spectrum of star # 5865 where the 
Mgl triplet is clearly visible. The spectra of two other peculiar 
hot stars are plotted for comparison. The spectra have been ob- 
tained for each star by summing the high-resolution ones used 
in Paper I, after correcting for RV shifts. 



fected to different extents and some are even not enhanced at 
all. The pattern of abundance ratios and its effect on the at- 
mospheric structure are still poorly understood. The approxi- 
mation of solar abundance ratio could introduce a systematic 
effect in parameter determination to an unknown extent, and in 
principle stars with different abundance patterns (due to differ- 
ent formation history) could beh ave differently in our a nalysis. 
Eitherway, there are indications (Edel mann et al.ll2006l) that, at 
least for EHB stars, the surface abundance pattern could almost 
be independent of stellar parameters, so any systematic effects, 
if relevant, should be the same in a wide range of tempera- 
tures a nd gravities . This m ight not be the case for cooler stars. 
In fact, IPace et alJ (120061) find a clear trend toward increasing 
heavy elements abundance with effective temperatures across 
the Grundhal jump (10000 < T eff < 14000 K). Therefore, 
at least in this temperature range, the deviations from solar 
abundance ratio (and then the differences with respect to the 
models) should be increasing with effective temperature. This 
could help for explaining our results, but the extent of the con- 
sequences on derived parameters is unknown, and in Fig. [7] a 
clear trend with temperature in this range is not visible. 

Our results for the three peculiar stars at 14000 K appear 
to be compatible with the hypothesis of helium-enrichment 
(as shown in Fig. 0J. However, in this case the low surface 
gr avity is directly lin ked to the higher luminos ity as shown 
by lCatelanl d2005h and lSweigart & Catelanl dl998l) . The masses 
derived from effective temperatures, surface gravities, and lu- 
minosities should therefore be well above 0.5 Mq. The low 
masses derived here suggest that the low gravities are at least 
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partly spurious. The same reasoning seems to exclude the hy- 
pothesis that these stars are an evolving or evolved object, 
as could be guessed looking at their position in the color- 
magnitude diagram. In fact, in this case again, the higher lu- 
minosities would compensate for the lower gravities, and the 
masses should agree with expectations. This is what we observe 
for star #18782 at U-Vx:-0.5, which shows both brighter lumi- 
nosity and lower gravity than other stars, but this is not pecu- 
liar for its calculated mass. Therefore we suppose that helium- 
enrichment or evolutionary effects can be considered good hy- 
pothesies that explain our observations for this star. 

It must be noted that the measured helium abundance can- 
not be properly disentangled from measured gravities because 
both have been measured with the same fitting procedure. 
Helium abundance s are mainly measured fro m He lines, but 
as investigated by Moehler & Sweigartl d2006l) . varying the he- 
lium abundance directly changes the width of Balmer lines, 
in the sense that higher helium abundances should produce 
broader lines. This rules out the possibility that low gravities 
are produced by an unnoticed helium enrichment, but indicates 
that underestimating gravity and overestimating helium abun- 
dance are two correlated effects. That these stars are systemat- 
ically brighter prevents concluding that the anomalous param- 
eter comes only from this weakness in the fitting procedure. 
These stars must be intrinsically different from fainter ones. 

Probably the best explanati on is that our result s indicate 
enhanced stellar winds. In fact Vink & Cassisil (12002b indicate 
that a non negligible wind should be expected for HB stars. 
They also show how it can change Balmer line profiles in such 
a way as to mimic a lower surface gravity. There is still a 
lack of detailed models that take the presence of wind into 
account, making it necessarily unaccounted for in our fitting 
procedures. Then, the derived low surface gravities could re- 
sult from this fault. The wind scenario is corroborated by the 
anomalous stars being brighter, so an enhanced stellar wind 
could be expected. Moreover, a wind could possibly counter- 
act the diffusion, which reduces the atmospheric helium abun- 
dance observed in subdwarf B stars to sub-solar values and thus 
justifies the higher helium abundances we measure. 

It is harder to find a good explanation for our results on 
the eight hot peculiar stars indicated as open squares in all 
plots. They were first spotted in Fig.|7]for their abnormally high 
masses, and then individuated in the color magnitude diagram 
as being fainter and/or redder than other stars observed in the 
same color-magnitude range. 

In this case we can assume that spectroscopically derived 
masses and photometric properties are independent. It is not 
generally true because magnitudes were used in the mass cal- 
culation, but as we saw, high masses cannot be ascribed to 
fainter magnitudes, and then the two effect cannot be related. 
The cause of peculiar masses must be other than photomet- 
ric data, and we would expect these stars to be randomly dis- 
tributed among the EHB population. Under this assumption, we 
performed a two-dimensional Kolmogorov-Smirnov (KS) test 
to explore the possibility that these stars are also photometri- 
cally peculiar. The two-dimensional version of the KS-test is 
n ot mathematicall y well-defined, but we adopted the algorithm 
of lPeacockl ( 11983b . which gives reliable results for all the cases 



of practical interest. We compared the distribution of the eight 
stars in the color-magnitude diagram with the whole population 
of stars in the range -1.25<(U-V)<-1 and 17<V<17.5. The test 
reveals that there is a 8.4% probability that the analyzed sam- 
ple is randomly drawn from the test distribution. This result is 
clearly not conclusive, but indicates that with high probability 
(91.6%) the eight hot peculiar stars indeed show photometric 
peculiarities and could belong to a sub-population of stars dif- 
ferent from the brighter and bluer bulk of EHB star population. 
It can be excluded that this result stems from a selection ef- 
fect in the target selection, because the same test, when applied 
to all 13 spectroscopically observed stars in this temperature 
range, gives a 25% probability of this sample being drawn from 
the test distribution. 

For the eight hot peculiar stars, the helium enrichment hy- 
pothesis is a possible explanation. In fact, in the T e ff-log g plot 
these stars agree better with polluted mode ls and, according 
to calculations by Sweigart & Catelanl l 1998 ). an increased he- 
lium abundance in the high temperature regime implies higher 
gravity. However, in this case their overly high masses again 
put this explanation in question, because higher gravities and 
lower magnitudes would balance out and lead to normal de- 
rived masses. 

As these stars are redder than their neighbors in the color- 
magnitude diagram, they might be binaries. In this case, con- 
tamination by the light of a cool companion could affect the 
parameter determination, and additional light would increase 
any derived mass. We therefore analyzed the spectra looking 
for signatures of the presence of cool companions, both in the 
low-resolution spectra used in this work and the higher reso- 
lution ones from Paper I, focusing on features typical of G-K 
type stars, such as the G-band and the Mglb triplet. We found 
no evidence of companions, except for star #5865, which shows 
a strong Mgl triplet, unusual for such a hot star (see Fig. 110b. 
Although this could indicate a binary system, it is just an iso- 
lated case, and the failure to individuate these features in all the 
other seven stars puts the hypothesis of binarity for these pecu- 
liar hot stars in doubt. In addition we did not find any indication 
in Paper I of radial velocity variations for any of these stars. 

In principle, we cannot exclude that their masses are ac- 
tually higher, but it appears to be quite improbable and hard 
to explain. It is simpler to think that the gravities are overesti- 
mated by fitting procedures, but even this conclusion has im- 
portant implications. As a consequence, we must deduce that 
models are unable to describe these objects and their spectra 
properly. Some difference between the model and the real star 
should be responsible for unreliably measured masses. But the 
same models give results in perfect agreement with theoretical 
expectations for "normal" stars, therefore the "peculiar" ones 
must be intrinsically different from the main population. This 
conclusion is strongly corroborated by the color-magnitude di- 
agram, where peculiar stars seem to be fainter and/or redder 
than the others. 

We point out that th e hot peculiar stars lie in concomitance 
of the Momany Jump (iMomany et al. I l2002h . a feature discov- 
ered first in NGC6752 and t hen also individuated in other 
clusters (Mo many et al. 12004b . The Jump is located around 
T e ff =23 000 K and is characterized by a sudden increase in lu- 
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minosity of the HB stars. As suggested by the author of the 
discovery paper, the presence of the jump could be due to deep 
changes in atmospheric structure at this critical temperature, 
with the onset of radiative levitation and changes in chemical 
abundances. We find no good explanation for our results in this 
direction, since we would expect brighter stars (and not fainter) 
to be deviating with respect to the models, but it is a coinci- 
dence that probably should be investigated further. 

6. Conclusions 

We measured atmospheric parameters (T e ff, logg, log ) 
and derived masses for 51 HB stars in NGC6752. Our re- 
sults agree in general with previous studies, but not always 
with theoretical expectations. Like previous works in literature, 
we find masses that are too low for cool (T e ff < 10000 K) 
and intermediate-temperature (12000<T eff < 15 000 K) stars. 
We find two groups of stars deviating with respect to the gen- 
eral behavior of the sample. The first group at 14,000 K shows 
masses that are too low, for which we consider a weak stellar 
wind (unaccounted for in the model spectra) as the most proba- 
ble explanation. The second group at about 27 000 K shows too 
high masses, for which we found no good explanation. These 
stars most probably (91.6%) do not belong to the general distri- 
bution of EHB stars in the color-magnitude diagram. We con- 
clude that atmospheric models, successful in reproducing the 
other EHB stars (for which we obtain masses in agreement with 
expectations), are inadequate for these peculiar ones. 
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